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Abstract 
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Dalitz analyses are introduced as the method for studying hadronic decays. An 
| accurate description of hadron final states is critical not only to an understanding 

of the strong coupling regime of QCD, but also to the precision extraction of CKM 
"P^ . matrix elements. The relation of such final state interactions to scattering processes is 

Qh| discussed. 

CD . 

1 Motivation 

><■ 

This serves as the introduction to following talks on Dalitz plot analyses. I will first remind 
you of the motivation for such analyses and discuss some aspects of how to perform them 
and I leave it to others to describe detailed results. 

We know that there is physics beyond the Standard Model but we do not yet know what 
this is. Insight is provided by precision measurements of the CKM matrix elements, and 
in particular the length of the sides and the angles of the CKM Unitarity triangle. Such 
measurements in heavy flavour decays probe the structure of the weak interaction at distance 
scales of O.Olfm. However, the detector centimetres away records mainly pions and kaons 
as the outcome of a typical process like B — > D(— > Kittc)K decay. Consequently, the very 
short distance interaction we wish to study is only seen through a fog of strong interactions, 
in which the quarks that are created exchange gluons, bind to form hadrons and then these 
hadrons interact for times a 100 times longer than the basic weak interaction. To uncover 
this basic interaction we need to understand the nature of these strong processes, or at least 
to model them very precisely. Indeed, the biggest uncertainties in determining the CKM 
angle 03, or 7, from the difference of D and D decays is due to our inability to model the 
final state interactions [H 121 E] . Dalitz plot analyses are the way to improve this. 
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What helps here is that the most common particles for 5's, D J s, and even J/ip's and (/>'s, 
to decay into are pions and kaons. These, being the lightest of all hadrons, have final state 
interactions that are common to all these processes. So each can teach us about the other 
and in turn heavy flavour decays are now the richest source of information about light quark 
dynamics. 
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Figure 1: Example of a Dalitz plot, here for D — ► K s tt + tt decay. 



As soon as one has more than two particles in the final state, the events are most readily 
pictured in a Dalitz plot, an example of which is shown in Fig. 1. This was invented by 
Richard Dalitz (known to everyone as Dick) more than 45 years ago [3]. Sadly Dick died 
earlier this year aged 81 and this talk is dedicated to him. For a three body decay, like 
D —> Ktttt, that we will concentrate on here, the data are plotted with the mass squared of 
the Ktt and tttt systems on the x,y axes, Fig. 1. Now the first thing to notice about such 
a plot is that the events are not uniformly distributed. The decay does not proceed by D 
decaying to the three body system Ktttt directly, rather there are a series of bands which 
show that the D likes to emit a K and then form a resonance like the p, which later decays 
into two pions, or the D emits a tt forming a if* (892), which then later decays to Ktt. Thus 
to describe the basic matrix element for the decay we need to know how to represent the 
vector mesons like p and K*. These, being reasonably long lived, have a magnitude and 
phase change across the resonance well described by a simple Breit-Wigner formula, with a 
pole in the complex energy plane on the nearby unphysical sheet. However, this represents 
only a small fraction of the events in the Dalitz plot. 
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The simplest way for any state of heavy flavour to lose mass is to emit a it or K and form 
a scalar meson. Having J PC = ++ this produces no change in angular momentum and so 
is almost cost free. Indeed, any system can lower its mass by forming a scalar with 1 = 0, 
since this has vacuum quantum numbers. This we can study in tttt scattering. There we 
see that the cross-section, Fig. 2, does not appear to have any structures looking like simple 
Breit-Wigners, rather it has a series of broad peaks with deep narrow dips between [Sj. The 
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Figure 2: Sketch of the modulus squared of the / = J = tttt — > tttt amplitude from Ref. 5. 
Resonances in the PDG tables [Hj are indicated. 

increase from threshold is related to the short-lived a and the dramatic dip at KK threshold 
to the /o(980). Neither can be represented by a simple Breit-Wigner and certainly naively 
summing Breit-Wigners would violate the conservation of probability. Consequently, we will 
need a better way to treat these broad and overlapping structures throughout the D and 
D s decay regions. Indeed data on these decays provide unique information about the strong 
coupling regime of QCD, which defines the Higgs sector of chiral symmetry breaking. Even 
if we do not care about this key aspect of the Standard Model, we nevertheless need an 
accurate way to parametrize its effects to determine the CKM matrix elements. 

We have seen that much of 3-body decays proceeds as a two stage process. This is 
embodied in the isobar picture, Fig. 3. There one assumes that parent particle P spits 
out a particle c, for instance, to produce a resonance R that subsequently decays into ab, 
particle c being regarded as a spectator as far as the final state is concerned. Such a picture 
immediately implies a modelling for the process ab — > ab, in which the same resonance R 
appears. One must check that these descriptions are consistent. The resonance must have 
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Figure 3: Decay of a parent particle P to abc. R is a resonance in the ab channel, which is 
the basis of the isobar model of sequential decays: P — > Rc, then R —* ab, plus contributions 
P — ► P'a, then P' — > 6c and P — > P"6 and R" — > ca with c, a and 6 as spectators, respectively. 

the same mass and total width, as well as partial width to the ab channel. This constraint is 
encoded in unitarity, which enforces the conservation of probability. Indeed, unitarity does 
not differentiate between the resonance and its background, but treats them both together. 
This is particularly important for short-lived states like the a, where such a distinction is 
totally semantic. Thus unitarity requires that if the particle c is a spectator then for instance 
the imaginary part of the amplitude for P — > 7itt(c) is equal to the sum shown in Fig. 4, 
where one sums over all kinematically allowed hadronic intermediate states. This means 
that if we represent the scattering amplitude for these intermediate processes ab — > h n by 
the matrix T, then the amplitude for the decay F is given by the vector equation 

F{P^ab{c)) = ]Ta„T(/i n ^ ab) , (1) 

n 

where the functions a n represent the basic coupling of P — > h n (c) for each hadronic inter- 
mediate state h n . 
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Figure 4: Unitarity for the ab system in the decay P — > (ab)c, where c is a spectator. The dashed 
lines denote the particles in the intermediate state are on mass shell. 

When c is a spectator these coupling functions are real. The hadronic amplitudes T can be 
conveniently represented by a fT-matrix form to ensure unitarity for the hadronic amplitudes 
(or by related forms I call the L- matrix if we are also to require the correct analyticity). At 
the lowest energies when ab is the only accessible channel this implies Watson's theorem jjj, 
namely that the phase of each decay amplitude (with definite isospin and spin), must equal 
that for the corresponding hadronic scattering amplitude for ab — > ab. 
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2 Semi-leptonic decays 



Now we can test this by considering the semileptonic decay D + — > K~7T + (fi + i> fl ), where 
clearly only the Kir system can have strong final state interactions. The dominant low 
energy signal is in the region of the if* (892), so FOCUS [S] looked at the mass region 
from 800 MeV to 1 GeV in two bins. If the vector if* was all there was then the angular 
distribution of the final state hadrons in the if tt rest system would be proportional to cos 2 9 
and so forward-backward symmetric. However a marked asymmetry is found, Fig. 5. This 
means the P-wave Kn amplitude must interfere with some other wave and at low energies 
that is the S-wave. It is here that many have argued lies the scalar k. What does experiment 
tell us? 
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Figure 5: Feynman graph of the semi-leptonic decay of a D + . Asymmetry in Kir angular distri- 
bution in the if* (892) mass region from the FOCUS experiment jHJ. 



Let us first look at Kir elastic scattering. This we learn about by studying high energy 
K 7r production with if -beams at small momentum transfers, where the reaction is controlled 
by one pion exchange. Then one can extract the cross-section for Kir — ► Kir scattering. This 
was done 25 years ago by LASS The cross-section shows clear if* (892) and if* (1430) 
peaks. By performing a partial wave analysis :9J, one confirms that these have spin-1 and 
2 respectively with the magnitude and phase of well-defined resonances. What we are in- 
terested in here is the .S-wave underneath these. This shows a broad rise in magnitude and 
phase from 825 MeV (where data begin) with a peak at 1400 MeV and phase rising through 
90°. This is characterised by the ifo(1430) with a width of ~ 300 MeV. There is no hint of 
a low mass k, which if it were simply describable by a Breit-Wigner, would have forced the 
phase to have already reached 90° by 800 MeV. What does D semileptonic decay tell us? 
FOCUS B\ found that their forward-backward asymmetry, Fig. 5, requires an S-wave phase 
of 45° in the if* (892) region, exactly as LASS has measured [0]. 

The statistics at FOCUS only allow a determination of the ifrr S-wave phase in the 
if* (892) region and then only in very wide bins. However, such D semileptonic decays with 
the event rates of CLEO-c and 5-factories should provide an accurate determination of the 
low energy if tt phase-shifts, in the same way as if e4 decays have done for irir phases. These, 
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when combined with dispersion relations and three channel crossing symmetry, have resulted 
in a rather precise determination of the position of the cx-pole JT0]- If we are going to be 
able to determine whether and where a k exists, then we will need precision information 
below the range accessed by LASS. Semileptonic D decays are the theoretically unambiguous 
way to go. 

3 Hadronic final states 

While we wait for that, let us turn to 3-hadron decays and let us focus on the Cabibbo- 
favoured D + —> K~7t + tt + channel. Let us consider this with increasing levels of sophisti- 
cation. First let us assume an isobar picture. The Dalitz plot is to be described by a sum 
of isobars in the three di-meson channels, Fig. 3. Since the 7T7t channel has / = 2, it has 
no known resonances and so is set to zero. The Kir channels are described simply by sum- 
ming Breit-Wigners with parameters taken from the PDG tables jH], the K*(892), ^(1410), 
-^2(1430) and of course the broad _K"q(1430). To this is added the 3-body interaction matrix 
element, which is presumed constant in both magnitude and phase across the Dalitz plot. 
The resultant fit to the E791 data [TT] is very poor and 90% of the decay is ascribed to the 
direct 3-body term. Since the Dalitz plot displays distinct 2-body structures, no wonder the 
fit is poor. The next step is to add a Breit-Wigner for another scalar, they call the k. It is 
just added. Its mass and width are determined by the fit [T2J to be M = 797 MeV, T = 410 
MeV. The fit dramatically improves. The direct 3-body fraction is down to a more believable 
14%, and the k contributes 43%. However modelling the scalar channel by a sum of the k 
and i^o(1430) Breit-Wigners is not only in violation of the conservation of probability, but 
implies a model for Kir elastic scattering in total disagreement with the LASS results |9, . 
We clearly have to do better. 

Rather than be tied to specific forms for the complicated S'-wave amplitude, Brian Mead- 
ows p"T|IT3] has tried something much more promising. With P and D-w&ve Kn interactions 
given as before by sums of Breit-Wigners from the PDG tables jB] , the S- wave is parametrized 
in terms of a magnitude and a phase in each bin across the Dalitz plot. Fitting gives the 
phases shown in Fig. 6. These are determined almost to threshold and, if they had greater 
precision and we applied the right tools for analytic continuation, might result in locating a 
k pole on the nearby unphysical sheet. So how do these phases compare with the Kir elastic 
S'-wave? 

In Fig. 6 the Kir phases from Brian Meadows' E791 fit ^3] arid the phases for K~tt + — > 
K~ir + scattering from LASS jU] are shown. The latter are absolute. Those from D-decay 
are relative, relative to the P-wave fixed to be 90° at 892 MeV. Consequently, we are free 
to raise the E791 phase up to be zero at threshold. We then see better agreement, with a 
common rising trend with mass. However, Fig. 6 shows these are clearly not the same. 
There are several possible reasons for this. K~tc + interactions involve both / = 1/2, 3/2 
components. In elastic scattering the relative strength of these is fixed by Clebsch-Gordan 
coefficients, while in D-decay this is determined by dynamics. We may have the prejudice 
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Figure 6: The points labelled LASS are the 5-wave phase for K~ir + — > K~ir + from analysis 
of the LASS experiment jS]. The solid data points are the K~n + S- wave phase from Meadows' 
analysis jlll I13j of the E791 results on D + — > i^ _ 7r + 7r + decay shifted from their original position, 
marked by crosses, by ~ 100°. 



that / = 1/2 is dominant, but that does not mean the 1 = 3/2 component is negligible. 
Thus if we assume Watson's theorem, which holds in the elastic region effectively up to Krj' 
threshold, one can determine the relative amounts of I = 1/2 and 1 = 3/2 contributions, as 
in Ref. 14. However, what I want to present here is something a little different. 




Figure 7: Unitarity for the Kir system in _D-decay in the elastic region. The dashed lines denote 
the particles in the intermediate state are on mass shell. 

So far we have considered the final state interaction of the K and a 7r where the second 
pion is merely a spectator. Let us now ask what happens if we try to include a subsequent 
interaction of the K with this pion. There has for long been a body of work on such 
multiparticle interactions particularly by Ascoli and collaborators 16 on 3 pion final 
states dating from the discovery of the a\ and its possible structures. Such multiparticle final 
states have been investigated by Anisovich et al. J7j in pp annihilations at LEAR. Much more 
recently Caprini ^H] has shown that one can deduce a unitarity relation with rescattering, 
implicit in these studies. This relation shows that for each partial wave amplitude the 
imaginary part of the D — > Kim is a sum of contributions, Fig. 7. 
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Figure 8: The solid curve denotes a fit to the LASS results in Fig. 6 consistent with Chiral 
Perturbation Theory |19j for the 1 = 1/2 S- wave Kir phase compared with the S- wave Kit phase 
from the Meadows' analysis of the E791 results on D — > Kirir decay. The upper line displays 

a preliminary calculation of how the elastic scattering phase incorporated in the first graph on the 
right hand side of the unitarity relation shown in Fig. 7 is modified by the rescattering corrections 
from the second (and third) graphs in Fig. 7. The curve closest to the E791 phases includes 
feedback. 



If there was just the first term in Fig. 7 then Watson's theorem would hold in the region 
of elastic unitarity. The second and third graphs in Fig. 7 give corrections. Since the Dalitz 
plot for D + — > K~7t + tt + is symmetric in the two Kit systems, and assuming the 7r + 7r + 
amplitude to be negigible, we can then use the Meadows' results to compute the rescattering 
corrections given by these graphs. There are some technicalities I won't go into here, but 
the crude result gives the phase of the S'-wave K~it + interaction in D-decay to be the upper 
line in Fig. 8. The corrected phase is computed at each datum of the Meadows' analysis 
and then these are simply connected by straight lines. The similar line a little bit lower 
(and closer to the data) involves including feedback to improve the phase variation across 
the Dalitz plot. 

These preliminary results are encouraging. However, details still need to be checked 
before they can be considered definitive. Nevertheless, this method holds out the prospect of 
using high statistics data on hadronic -D-decays from BaBar and Belle to determine the near 
threshold Kir phase-shift with precision, independently of summing higher orders in Chiral 
Perturbation Theory. So while there is no k(900) (2H|, there may be a scalar much closer 
to threshold deep in the complex plane. From the Roy equation analysis of tcti scattering 
we now know the position of the cx-pole rather precisely. All analyses should find the same 
result. The variation between treatments discussed in Refs. 20, 21 is unacceptable in the 
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era of precision physics. The application of the rescattering corrections to J/ip — > urna may 
well, with higher statistics accessible at BESIII, bring the mass and width from a simple 
Breit-Wigner treatment [22] i n line with the true cr-pole position ^Uj. That is the challenge 
for BES. 

A precise description of these hadronic final state interactions, especially those in scalar 
channels, is essential to reducing the uncertainty in the CKM triangle. This reduction 
is crucial to learning about physics both within and beyond the Standard Model. Dalitz 
analyses are at the heart of this programme, as others will describe. 
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